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We show that statistics is crucial for the instability problem derived from higher time derivatives. 
In fact, and contrary to previous statements, we check that when dealing with Fermi systems, the 
Hamiltonian is well bounded and the quantum states are normalizable. Although, ghost states 
are still present, they do not affect unitarity under certain conditions. We first analyze a quantum 
oscillator involving Grassman variables and then we generalize it to a Dirac field. Finally, we discuss 
some physical implications. 
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In recent years there has been growing interest in 
field theories with higher derivatives. Nonlocal effects 
in string theory [Jl , quantum gravity and noncommu- 
tative geometry |3l| provide a natural landscape for the 
incorporation of higher derivative terms. Also, attempts 
to incorporate such terms have been given in the con- 
text of dark energy Lorentz violation radiative 
corrections [ij and regularization 7] . 

However, according to current belief, field theories with 
higher time derivatives lead to an unavoidable energy in- 
stability. In general, this statement is a consequence of 
Ostrogradsky theorem A paradigmatic example is 
given by the well known Pais-Uhlenbeck model Q where 
it is possible to show explicitly that the energy is not 
bounded below. One way to avoid Ostrogradsky theo- 
rem consists to consider a non normalizable ground state. 
However, this alternative not only has the drawbacks of 
containing this non normalizable CTOund state, but also 
yields states with negative norm [3, [13] ■ Although, 
those ghost states may cause a loss of unitarity, under 
certain conditions, they do not imply a serious problem 
in order to define a sensible theory 111!]. However, either 
the instability or the non normalizability invalids any at- 
tempts of a perturbative expansion. A third approach 
to deal with this problem has been made by redefining 



the inner product by using VT symmetry [12l |. Others 
include perturbative iterations [13] and performing com- 
plex transformations [l4i] . 

So far, nobody has explored the role played by statis- 
tics on this issue. In this work, we propose a model 
somehow equivalent to the Pais-Uhlenbeck model, but 
using fermion instead of boson variables. We show that 
our model mimics some features of the Pais-Uhlenbeck 
model: one frequency has no perturbative behavior, the 
phase space is increased, unitarity may be lost when an 
interaction is turned on. However, the instability prob- 
lem does not show up. The mechanism that allows this is 
essentially the same as the one that stabilizes the Dirac 
sea in the electron theory. 

To begin with, let us recall the Pais-Uhlenbeck La- 



grangian 
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which is basically the harmonic oscillator with frequency 
UJ plus a higher derivative term. The Hamiltonian of the 
system after quantization can be written as 



Hpu = k- 



k.b^b- 
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with fc^ = 1=F^/1 + Ag) positive frequencies depend- 
ing on g and lo and a, a) , b, the standard creation and 
annihilation operators. The second term produces arbi- 
trary negative energy states as can be seen by acting W 
on the empty wave function (defined by a$o = ^'I'o = 0) 



^0 ^ N exp 



'2ik, 
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The fact of having energies not bounded below, invalids 
the perturbation formalism. Indeed, the expectation val- 
ues of some observables depending on q, blows up when g 
goes to zero. An alternative proposal would be to redefine 
the vacuum a^g = ^^'^'o — 0' which solves the stability 
problem. However, this new vacuum state is not normal- 
izable, and negative normed states emerge. Our proposal 
is to change commutators with anticommutators, namely, 
to work with fermions instead of bosons. 
Let us then consider the fermionic model 



Lp — ^—ipip + iipip — utipip. 

UJ 



(4) 



As before, this is the (fermionic) harmonic oscillator 
{g — 0) plus a higher derivative term. Without loss 
of generality we can assume that uj and g are positive 
constants. The equations of motion are. 



tp = i—ip tp, 
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and its conjugate. If we try a solution like ■0 = % e 
we find that, A — A— — — = 0, whose solutions can be 



expressed in terms of positive quantities a;±, 



u}± = ±X± = cj, 
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where a;+ > w_ > 0. 

The Hamiltonian and the commutation relations which 
reproduce the equations of motion are, 



H = —iptjj + ujipip-> 
ui 

and the non vanishing anticommutators 



(7) 



{^,V} = «-, {i^J}^-"^. (8) 
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Using the Schroedinger representation in terms of Grass- 
man variables and their derivatives, we have, 

-;- uj d UJ - 
5 dip 2g 



(jj d ijj 
5 o^} 2g 
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And the Hamiltonian in this representation is. 



H = AX 

g dip dip 



uj[l + ^]ipip 



2g {"^diP "^diPj 2g 



(10) 



The first and second lines of this equation commute to 
each other, so it is very easy to find the eigenf unctions. 



V2 
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Note that for g > they are ordered in increasing energy, 
and of course, the system is bounded. So, the vacuum will 
correspond to the lowest energy state, i.e., the state $o- 
However, the shift of energy respect to the lowest energy 
state is what physically matters. After this shifting we 
obtain AEq = 0, AEi = w_, AE2 = w+ and AE3 = 
w+ + w_ which are all positive. Thus, we note that the 
energies of the first two states go to and to w as 5 
goes to 0, respectively. However, the last two energies 
are infinitely large as g tends to zero. 

Let us see what happens with the normalization. If we 
define the Berezin measure as, / dipdipipip — +1 we can 
see that, ($o|*I'o) = ('I'll'I'i) = 1 which are normalizable 
states (if g > -1/4). However, ($2! $2) = ('I'sl'&a) = 



— 1 are states with negative norm. It is worth noting 
that, since the system is bounded below, in the limit of 
5—^0, the ghosts states live in an arbitrary high energy 
level, and we should expect this sector decouples from 
the theory. 

Even in the case when g is finite, unitarityis pre- 
served. To see this, let us define the new inner product, 
where g (a|/3) = {a\G\(3), where G = exp(i7riV+), being 
iV+ — c\_c^ the number of ghosts operator. This new 
inner product is the usual positive defined one. This is 
to say, we can write {a\(3) = {a\G\(3), and then, the time 
evolution of amplitudes are. 



idt{a\f3) = {a\{GH - H^G)\P) 
= {a\{H -GH'^G)\13) 



(12) 



where W is the usual hermitian conjugate of H . Then, 
in order to obtain probabilities conserved we need H = 
GH^G. In the literature, ^15|], this is called pseudoher- 
micity condition, and, under certain circunstances, it is 
possible to obtain a well defined unitary theory [3| ■ This 
is the case for our modal, where G^ — I, W = H and 
[G, H] ^ 0. Then, we do not lose unitarity. Adding in- 
teractions can change the situation, but there are still 
chances to keep unitarity. In any case, when g is small 
enough, so that the ghost sector can be considered as 
decouple from the theory, we can perform pertubative 
calculations in a safe way. This last fact, is contrary to 
the bosonic case, where either the instability problem or 
the non normalizability of ground state invalids any per- 
tubative formalism. 

In order to clarify a little more what is going on, and 
to make a closer contact to the QFT formalism, let us 
consider the new variables. 



= a{iip + uj+ip), ~ a{—iip + to+ip), (13) 



and, 



c+ — a{—iip — uj-ip), c\_ = a{iip — uj-ip), (14) 



where a — ^ . It is easy to show that the non 

vanishing anticommutators are, {c_,c^} = 1, {c+,c^} = 
— 1. Then, the operators c„, cl^, c+ and arc standard 
creation and annihilation operators, but c+ and create 
and annihilate states of negative norm. 

The original variables can be written in terms of these 
operators as 

J 



4, 



VTT4^ 

so the Hamiltonian turns out to be 

Hp = uj-c^_c- — w-|_ct c+ — 



(15) 



(16) 
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This expression can be compared with the Pais- 
Uhlenbeck Hamiltonian [2] They are the same expres- 
sions, except for an irrelevant constant. However, as we 
see, the spectrum is bounded below, and the vacuum is 
normalizable. In fact, the true vacuum is the state anni- 
hilated by c_ and c+, namely. 



where the index s stands for the spin quantum number. 



-*o^o^(^-^^).o^0, (1.) 

which agrees with the result above. Morever, the rest of 
the states are $i = cL'&o, ^2 = c^^^o and $3 = cLc+$o 
with their corresponding energies. Then, the operator 
is creating ghosts states in the system, even though the 
problem is bounded and hence stable. 

This result can be generalized for terms with higher 
time derivatives whenever the frequencies hold real. Fur- 
thermore, the creation and annihilation operator formal- 
ism we have just seen, allows to generalize straightfor- 
wardly our quantum mechanical model to a fermionic 
QFT system with higher time derivatives. 

In the remaining of this work, we deal with this gen- 
eralization to a field theory framework. The obvious co- 
variant QFT generalization of our model is given by the 
Lagrangian density 



m)^ — —ip □ tp, 



(18) 



where 5 is a dimensionless positive coupling constant and 
A is an ultraviolet energy scale. This scale is the finger- 
prints of some fundamental theory and hence our model 
must be considered as an effective theory valid only for 
energies far below that scale. 

The equations of motion followed from (|18p are 



(i^ - TO - □ )?/; = 0. 



(19) 



The solutions of the equation are given in terms of plane 
waves ip{x) ~ J d'^pe~'^^'^w{p) where w{p) is a Dirac 
spinor satisfying 



9 2 



'w{p) ~ 0. 



(20) 



Nontrivial solution of this equation imply p^ = irij. 
with, 



±1 + ^1 + 4^ ^ 

TO± = ^A. 
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It turns out that the w spinors are the free Dirac solutions 
with masses m±, namely 



m± u±(p±) = 0, 
hTO± 'y±(p±) = 0, 



(22) 



and p± = {E± , p) with E± — ^p'^ + m\ . As expected, 
we have in the limit g — > one regular solution rri- ^ m 
and the other to-|_ going to infinity. Also, it is worth 
noting here that, even in the case when to is zero, we 
still have a non vanishing to-(- — A/g. 

Following the same spirit as in the last part of our 
quantum mechanical model, we can define the fields. 



■0_ = I3{i0 + m^)ip, 



(23) 



with 13 — ^ rn++m ) " Lagrangian density in terms 
of these fields is 

£ = ■0_(i^ - TO_) V'- - '0+(j^ + m+) (24) 

Now it is easy to see that the Hamiltonian is 

H — J d^x ■ J + m^ )ip 

- V;+(zV-7-TO+)7^+] , (25) 
with the non vanishing anticommutators, 
{V'-(x),^L(y)} = -{V^+(x),vA(y)} = <5^(x- y). (26) 

Now, decomposing these fields in terms of plane wave 
solutions we find that 



^4 



d^p 



^ J (27r)3 

[bi {p)e-'P- -^ui (p) + {p)e'P- -""vl (p) 
X ^ f d^p TO+ 

[b'_^{p)e-^P+--uXip) + d;(p)e'^f+-<(p)" 



(27) 



(28) 



Notice the change respect to the standard decomposition 
in ?/'+, swapping the b's and d^'s operators by b^'s and 
d's. The non vanishing anticommutators are the standard 
ones for 6_ and d_, that is, 



{51 (p), felt (k)} ^ 5-^3(p_k), 
K-(p),d':t(k)} - 5-53(p_k), 

but for the 6+ and (i+ we have, 

{6;(p),6;^(k)} ^ -5-53(p_k) 
{d;(p),d;t(k)} = -<5-<53(p-k) 



(29) 



(30) 



The minus sign is revealing that these modes are ghosts, 
i.e., these operators create and destroy negative norm 
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states. The Hamiltonian in terms of these operators is, 
except for an irrelevant constant, 

H ^Y. f rf'p(^^-(&l^(p)&l(p) + rf-^(pK-(p)) 

s 

+E+ (6^ (p)fe; (p) + (p)d; (p))) . (31) 



The vacuum is the state annihilated by all operators 5j_ 
and dj-. This state is the true lowest energy state, and it 
is normalizable and well defined. Hence, again we have 
a stable theory, but with negative norm states above en- 
ergies of order greater than m+ ~ A. As we discussed 
above, this fact ensures that if we add any known inter- 
action, the theory would still be well defined, as long as 
the energy scale involved in the physical processes are 
lower compared to the ghost mass scale. 

As it was discussed for the QM model, our QFT the- 
ory will not have unitarity problems. Only by adding 
interactions they can come up. However, the loss of uni- 
tarity due to the ghost sector, will be negligible as we are 
much below the ultraviolet scale. Even more, the stan- 
dard unitary theory is recovered when g tends to zero. 
This situation is similar to the one pointed out in 



but in our case we have a stable vacuum and we can make 
sense of the theory even without going to the Euclidean 
space. 

Finally, to calculate propagators, it is useful to have in 
mind that the original fields are combinations of tp- and 
namely. 



And the propagator is, 
1 



-[i^.{x)-iP+ix)] (32) 



S{x) 



(S^ix) - S+ix)) 



(33) 



where S± are the standard propagators for Dirac parti- 
cles with masses m±. The minus sign in front of 5+ is 
another signal of the ghost sector given by the -I- modes. 
This propagator resembles the Pauli-Villars regulariza- 
tion, with to+ as the regulator which render ultraviolet 
divergent integrals finite. However, this regulator works 
differently from the Pauli-Villars formalism in some sit- 
uations. For instance, it produces new vertices when we 
introduce interactions when we keep gauge invariance. 

Summarizing we have found that Fermi statistics can 
fix the problem of stability in higher time derivatives sys- 
tems, giving a counterexample of the Ostrogradsky's the- 
orem. However, negative norm states are still present. 
These ghost sector can be controlled either by restrict- 
ing ourselves to the physical sector or by keeping them 
far abouve the UV scale. In both cases, unitarity is pre- 
served as long as the relevant physical processes energy 



scale is far below the UV scale. And hence, it is well 
justified to add local interactions and perform perturba- 
tive calculations in our theory. Somehow, the situation is 
similar to the standard QED, where perturbation treat- 
ment are justified only for energy scales far below the 
Landau pole. Just above this pole, the fine structure 
constant becomes negative and amplitudes like 'vacuum 
to vacuum' might be negative, leading to possible ghost 
states. That is a signal of the fact that QED can be 
considered as an effective theory valid only for certain 
energy sclaes. It would be interesting to study physical 
consequences derived from these new terms in relevant 
theories like QED-like or symmetry broken theories with 
fermions 1171. 
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